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Abstract. We present both theoretical and experimental results on fluorescence 
of single defect centers in diamond nanocrystals embedded in a planar dielectric 
microcavity. From a theoretical point of view, we show that the overall fluorescence 
collection efficiency using moderate numerical aperture microscope objective can be 
enhanced by using a low quality factor microcavity. This could be used in particular 
for low temperature applications where the numerical aperture of collection microscope 
objectives is limited due to the experimental constraints. We experimentally 
investigate the control of the fluorescence spectrum of the emitted light from a single 
center. We show the simultaneous narrowing of the room temperature broadband 
emission spectrum and the increase of the fluorescence spectral density. 



PACS numbers: 78.67.Bf, 42.50.Ar, 42.55.Sa, 78.55.Hx 
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1. Introduction 

Single nitrogen-vacancy (NV) defect centers in diamond is a promising candidate 
for solid state quantum computing and quantum information processing (QIP) 
applications [lj. It has been shown that NV center can be used as a reliable source 
of single photons [2] for quantum key distribution (QKD) [31 H] and single photon 
interferences [5]. More recently, its potentialities as a spin qubit with a very long spin 
coherence time [Q] have been successfully demonstrated in the context of the coherent 
coupling of NV - electronic spin to single neighbouring nuclear spins [7J [8]. It takes 
advantage of both the spin coupling to optical transitions and the optical readout of the 
associated electron spin resonance [9]. The spin and optical properties of a single NV 
center were also used to probe magnetic fields at the nanoscale [TUJ CLE] • 

In both QIP or QKD applications, the efficiency of the collection of the radiated 
photons is crucial. The rate of the qubit exchange in QKD or the efficiency of 
optical read-out in QIP strongly depends on the external quantum yield. To control 
their spontaneous emission rate or to increase their fluorescence collection, solid state 
individual emitters can be embedded in an optical microcavity (T2J H5] . For NV 
center systems, different configurations of coupling with a microcavity have been 
investigated. Some efforts have been made in designing diamond photonic crystals 
microcavity [HI [151 HE]- Microcavity with a quality factors Q around 600 have been 
fabricated in a thin photonic crystals membrane [I7J. Due to the difficulty to obtain 
monolithic microcavities from bulk or nano crystalline diamond [IE] others approaches 
based on hybrid structures have been used. Normal mode splitting has been observed 
in a system consisting of diamond nanocrystals coupled to the whispering gallery mode 
(WGM) of a silica microsphere [19J . The coupling of a single NV center to a high-Q factor 
WGM resonator has been observed in silica microdisk [20] or microsphere [211 122] • By 
using a fiber taper, a single diamond nanocrystal has been positioned onto a silica WGM 
toroid which demonstrates the possibility to manipulate non-classical emitters in order 
to couple them with high-Q microresonators [23]. GaP WGM microdisk have also been 
used to efficiently collect at the chip scale the fluorescence of NV centers [21] . To reach 
very low mode volume, an optimization scheme of Ti02/Si02 micropillar bragg cavity 
has been proposed [25]. Alternatively, hybrid structure involving photonic crystals in 
Silicon Nitride [26J or GaP [27] enabling coupling to NV centers have been proposed 
to reach high Q-factors. In this context the controlled coupling of a single diamond 
nanocrystal to a planar Si3N 4 photonic crystal heterostructure [28] and the emission 
spectrum modification of a single NV center embedded in an opal photonic crystal 
structure made of polystyrene microsphere [29] were experimentally demonstrated. 
Finally one-dimensional photonic nanostructures consisting of diamond nanowires can 
also be used to enhance the single photon emission from NV centers [30] . 

In most of these previous cases, the goals were to obtain simultaneously a high- 
Q factors and a small mode volume in order to enhance the Purcell factor [31] and 
reach the regime in which the emitter spontaneous emission rate is strongly affected 
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by its coupling to a cavity mode. This results in enhanced spectral properties of the 
single-photon source well adapted to GHz repetition rate QKD |32j. 

However, the emission rate can be not only strongly modified but also the radiation 
pattern of the emitter [331 131] • The thorough investigation of the modification of the 
fluorescence spectrum and the radiative pattern by photonic structures is more recent, 
in particular for broadband emission. In that later case, the dispersion associated to 
its dielectric environment leads to a renormalization of the emitter spectrum, which is 
such that a simple planar microcavity can strongly modify the emission spectrum of an 
organic dye [351 EH1 Ell [38] or of a colloidal CdSe/ZnS nanocrystal [39J. 

In this work we show that the use of a planar microcavity lead to an enhancement 
of the collection efficiency of the single NV center photoluminescence in practical 
configuration of QKD or QIP. We also demonstrate that the microcavity can be used 
to tune the emission spectrum. The article consists of two theoretical and experimental 
parts. The first part begins with the description of the planar microcavity structure. We 
then present the calculation method used to determine the fluorescence power collected 
by the microscope objective from a single NV color center located either on top of 
a high-reflectivity mirror or inside the microcavity. We consider both cases of the 
monochromatic emission (standing for cryogenic operation of the NV center emission) 
and the room-temperature broadband emission. In the second part, we describe the 
experimental setup and the observations showing the modifications of the emission 
associated to the coupling of a single NV color center in a nanodiamond to a single 
longitudinal-mode of the planar dielectric microcavity. 

2. Calculation of the total fluorescence power P tot collected from a single 
emitting dipole coupled to a planar Fabry-Perot microcavity 

The collection of light from a single emitter embedded in a dielectric material with a 
high index of refraction, as NV color center in bulk diamond or semiconductor quantum 
dots, shares similar issues than those of semiconductor light emitting diodes (LED). 
They mainly come from the light trapping due to total internal reflections. A lot of 
theoretical and experimental studies have shown that light collection efficiency of LED 
can be enhanced by using a planar microcavity [101 HI] . This concept has been adapted 
to high index of refraction solid-state single-photon sources [13] . 

In the case of NV color center in diamond nanocrystals it has been pointed out that 
refraction is irrelevant due to the small dimensions of the nanocrystals compared to the 
emission wavelength and that the emitter can be assumed as a point source emitting in 
the surrounding dielectric medium [12]. This feature enhances the collection efficiency. 

In this section we will show that the use of a resonant microcavity combined with 
the nanocrystal approach can also increase the collection efficiency for NV defect centers. 
By using a high reflectivity rear mirror the isotropic radiation is redirected in the same 
half space [IB], while the microcavity simultaneously narrows the emission spectrum |44j, 
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both effects resulting in a larger spectral density of the collected light intensity . Finally 
for an appropriate choice of the emitting dipole position and orientation, one can also 
benefit from the local electromagnetic field enhancement [40J. 



2.1. Microcavity - single- emitter structure 

The single-photon source described in figure [U consists of a nanocrystal embedded in a 
thin polymer layer deposited on a dielectric high reflectivity (HR) Bragg mirror made 
of 16 periods of two materials of high and low indices of refraction |45j . 
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Figure 1. (a) Scheme of the single-photon source microcavity consisting of a diamond 
nanocrystal containing a single NV center, positioned in the air gap between the 
two mirrors of a planar Fabry-Perot microcavity. The lower high reflectivity (HR) 
Bragg mirror is composed of a fixed number of 16 periods, whereas the upper low 
reflectivity (LR) mirror is made of N periods. Referring to the figure annotations, 
the thicknesses values of the different layers are: h& = 134 nm, Ia = 171.7 nm, 
hi = hi = /ib = 73.1 nm, h = I2 = 113.7 nm, he — 75.3 nm, = 110.1 nm and 
/d = 117.1 nm. The single NV color center is assumed to be located on the top of 
the polymer layer of thickness zn — 25 nm. do is the adjustable air gap thickness. 
The low index of refraction (ri£ — 1.4745 at A = 680 nm) layers are marked L and 
the high (m, = 2.2925) index layers H. The refractive indices of the polymer and of 
the substrate are respectively n p = 1.49 and n s = 1.4556. We also give the optical 
intensity distribution |-E(z)| at A^y x for a plave wave excitation in normal incidence. 
The calculation have been made using the transfer-matrix method [46] for (a) the HR 
mirror and (b) the microcavity with An = A™y x an( l N = 4. 
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The cavity is obtained using a second low reflectivity (LR) mirror made of N < 16 
periods of the same materials. The resonant wavelength of the microcavity can be 
tuned by changing the air gap do between the LR Bragg mirror and the HR mirror. The 
thicknesses of the different layers of the Bragg mirrors have been chosen to obtain a 
central wavelength Ao around A™y x = 680 nm corresponding to the peak emission of the 
NV~ center at room temperature. The structure has also been optimized to maximize 
the electric field amplitude in the polymer layer containing the nanocrystals. We give 
in figure [U the electric field amplitude distribution for a normal incidence plane wave in 
the two cases we extensively study in this paper: (a) HR mirror without second mirror 
(b) microcavity with iV = 4. 
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Figure 2. Graphical representation of the Bragg mirror and microcavity dispersion 
relations (left) and of the transmission spectrum (right). TE (resp. TM) polarization 
case is displayed on the negative (resp. positive) fen side, a = h + hi is the period 
of the Bragg mirror. The propagation is forbidden in the white areas and allowed in 
the pink areas. The planar microcavity (MC) dispersion curves are also represented 
(in blue) assuming an interference order m c = 3 for Ao in normal incidence. We have 
added two light lines: the one associated to a full half-space collection (numerical 
aperture N.A.=1, red dot line) and the other for a restricted collection of light with 
N. A. =0.6 optics (grey line). Right graph: transmission spectrum T(u>, 9, TM) of the 
microcavity (N — 4 case) in dB units at either normal incidence (9 = 0°, blue curve) 
or in the case of inclined incidence with sin# = 0.6 and TM polarization (red curve). 
These curves display dips at their resonant wavelengths. The dashed region on the 
left graph indicates the range of frequencies allowed in the case of N.A.=1. Note that 
for N.A.=0.6, only the frequencies shorter than the resonance frequency are allowed, 
which explains why the experimental spectra displayed on figure [8] and figure [9] are 
asymmetric, with a steep falling edge on the high frequency side. 

Figure |2] represents the dispersion relation u)(k\\) of the Bragg mirrors where u is 
the angular frequency and k\\ is the in-plane component of the wavector k jl6]. Due 
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to the field penetration in the dielectric mirror the minimal interference order for the 
planar microcavity is given by m c = \n/An\ + 1 with An = — rig, n = n^ni/n^ and 
n & = (neh + n^hi) / (li + hi) where |_^J is the floor function of x [10]. With our refractive 
index values we obtain m c = 3 which is used to obtain the microcavity dispersion curves 
displayed in figure [2] as solid lines. 

The dispersion relation is obtained noting first that k = ^Jk 2 + k 2 z . Then, if we call 
L cav and n cav the effective length and refractive index of the cavity taking into account 
the penetration depth in the Bragg mirrors, we obtain 2k z L CSiV = 2m c ii for a resonant 
mode and thus: 

c /, 9 / m c 7T 

' 'cav y \ -^cav 

Neglecting the polymer layer effect, the two quantities L cav and n cav can be estimated 
by solving the resonance equation (for d w 135 nm): 

m c y = 1x4 + n a (L ca , v - d ) = n cav L cav . (2) 

We notice that with the microcavity design, the dispersion curve referring to m c = 3 
is fully included in the first forbidden band of the Bragg mirror for both TE and TM 
polarizations and for the maximal numerical aperture considered here (i.e. N.A.=1, 
corresponding to the light line) and thus also for the value N.A.=0.6 used in the 
experiments. This feature is important for room temperature experiments since in 
this case the fluorescence spectrum of NV centers is broad (full width at half maximum 
FWHM?^ 75 nm) and therefore requires a broadband operation range. We also show the 
transmission spectrum T(u, 8, TM ) obtained using the standard transfer matrix method 
[46J for a microcavity with N = 4 both for the normal incidence, and for sin 9 = 0.6 TM 
polarization cases. Note that due to the unbalanced mirror reflectivity, the resonant 
transmission is weak even though we assumed lossless materials. 

This simple linear analysis [13] allows two properties of the microcavity to be 
revealed: i) for the numerical aperture N.A.=0.6 used in our experiment it is possible 
to efficiently collect light coupled to the resonant mode of the cavity, ii) in the case 
of a broadband emission, the blue-shift of the cavity resonance explains the spectrum 
asymmetry experimentally observed (see Section [3]). 

2.2. Calculation method of the total collected power P tot radiated by a single NV color 
center inside the microcavity 

It is well established since the early studies of NV~ color center [27], that it has two 
orthogonal incoherent electric dipole transitions oriented in the plane perpendicular to 
the [111] crystallographic direction (along the N-V axis), as exemplified in the context 
of recent studies on polarization selective excitation of this color center [IHl H9] . 

In the following we will model the NV~ center by two identical orthogonal 
dipoles (e, ej_). The dipoles are assumed to be located in O. e orientation is defined 
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by two angles G and $ relative to the Bragg mirror plane (x,y), such as e = 
(sin cos $, sin sin <£>, cos 0) in the sample plane coordinate system as represented 
in figure El Thus the second dipole can be written: = Vi cos a + v 2 sin a with 
V! = (— sin $, cos $, 0) and v 2 = (— cos cos — cos sin <£>, sin 0). The most 
favorable configuration to efficiently collect the photoluminescence of the NV~ color 
center corresponds to (0 = n/2, a = 0) whereas the worst configuration is obtained for 
= 0. In the experiment we selected the brightest emitters, which are the one having 
the strongest in-sample plane dipole components. Therefore, in the following we restrict 
our calculation of the collected power P tot to this configuration (0 = ir/2,a = 0). To 
calculate the power radiated by the two dipoles in the upper half space (z > 0), we used 
the transfer-matrix method of reference |50j . 




Figure 3. (a) O is the position of the orthogonal (e, ej_) dipoles that we associate 
to the NV~ color center transition. The orientation of the component e is described 
by two angles (0, $) relative to the Bragg mirrors plane (x,y). k + is the wavevector 
of the field in the upper half space emitted in the direction (9, <j>). (b) a refers to the 
orientation of ej_ in the plane (vi, V2) orthogonal to e. (c) k_ is the wavevector of the 
field emitted in the same direction but in the lower half space. The unitary vectors 
are noted s± (resp. p±) for TE (resp. TM) polarization. 



For a given direction of radiation characterized by the wavevector k + , with 
coordinate k + = /c(sin# cos0, sin# sin 0, cos#), the wave can be either TE-polarized 
along s + or TM-polarized along p + . s + is the unit vector orthogonal to the incidence 
plane (u,z), whereas p + is calculated by: 

P+ = x k + . (3) 

The wavevector of the field related to k + propagating in the opposite direction 
is noted k_ = A; (sin 9 cos <fi, sin 6* sin 0, — cos^). We also have s_ = s + and the same 
relation as Eq. ()3]) holds for p_. According to the method given in reference [50J and 
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using the dipole description of references [3~J] and [5T] the source terms for the two 
propagation directions are: 

A s + = e ■ s + , A\ = e • p + , A s _ = -e ■ s_ and A p _ = e ■ p . (4) 

We have the same relations for e_i_. Using these expressions for the optical sources 
it is possible to calculate the spectral density of the radiated power per solid angle 
d 3 P/(dnd\) for a given k+ by summing the power emitted by the two incoherent dipoles 
[50] . Following the method described in reference [37], the full radiated power collected 
by the air microscope objective of numerical aperture N.A.=sin 6u in a spectral window 
of full width AA is then evaluated by : 

Ptot= f(X)dX / #/ S m6d6 (5) 

Jx -A\/2 JO JO CL\laA 

where /(A) is the free space emission normalized spectrum function which will be 
assumed for simplicity to have a Gaussian profile centered at A^y x with a full width 
at half maximum B. 



2.3. Calculation of Ptot 

Our aim is to show that the light extraction from a single emitter can be increased 
using a low finesse microcavity in practical cases. Consequently we will take as the 
reference the case of dipoles located on a single mirror. The advantage of this last 
configuration with respect to free space is twofold: the mirror redirects the fluorescence 
in the same half space and if the dipoles are located in a field antinode, their emission 
rate is increased [401113]. In the following we consider two cases: 1) the monochromatic 
emission of NV~ at its zero-phonon line (ZPL) wavelength Azpl = 637 nm, which 
corresponds to cryogenic operation, 2) the broadband emission case, associated to room- 
temperature photoluminescence. 



2.3.1. Monochromatic case First we focus on the low temperature configuration. In 
this particular case, the fluorescence spectrum of the emitter is assumed to be quasi- 
monochromatic and the ZPL has a bandwidth around B = 0.3 nm j52j [53]. For 
simplicity and because our cavity has been optimized for a resonant wavelength 680 nm, 
we have used for the calculation A^v* as the emission wavelength instead of Azpl- 

In order to collect the light across the fused silica substrate of the LR mirror in the 
case of the microcavity configuration, we use a long working distance objective with a 
numerical aperture N.A.=0.6. The reference configuration consists of a HR mirror and 
a microscope objective with N. A. =0.85. We used this low value of N.A. because it is 
suited to cryogenic setup [52] • The simulation results are shown in figure HI In the case 
of low finesse cavity (N = 1) and a detuning 5 = Xq — A^y* = 30 nm, we observe an 
increase of about 50 % in the collection efficiency with respect to the reference case. 
This can be increased by using a better finesse (J 7 ~ 60) cavity (N = 4). Note that 
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Figure 4. Full radiated power P tot in the case of a narrow bandwidth emitter (single 
NV~ color center emitting at A™y x ) as a function of the numerical aperture of the 
collection objective. The calculations have been carried out for a single HR mirror and 
in the case of the microcavity for two values of N and different values of the microcavity 
mode - emitter wavelength detuning 5. The results are normalized to the value of Ptot 
corresponding to the HR mirror case and at the maximal numerical aperture collection 
(N.A.=0.85). 

in this case, the cavity must be slightly positively detuned (5 = 10 nm) to obtain the 
maximal collection increase of about 80 %. This positive detuning allows the spectral 
overlap between the emission spectrum and the microcavity dispersion to be optimized 
as illustrated by the hatched area in figure EJ For 5 = we do not observe any increase 
in the extracted fluorescence for the maximal numerical aperture. 

2.3.2. Broadband emission At room temperature, the coupling between the radiative 
transitions of the NV~ color center and the diamond matrix phonon broadens the 
emission spectrum yielding a FWHM linewidth of about B = 75 nm due to the phonon 
replica. 

Figure 13(a) gives the calculation of the full radiated power with no additional 
filter whereas figure El (b) represents the same calculations with a AA = 40 nm band- 
pass filter whose resonant wavelength is chosen in order to increase the overall radiated 
power for the maximal numerical aperture. This filter could be useful for example 
in open space QKD applications where it will reject a large amount of background 
light [4J. As it was the case for monochromatic emission the enhancement is better 
for N = 4 and for a positive detuning. The maximal enhancement is weaker than 
in the previous configuration since in room temperature experiments it is possible to 
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Figure 5. Full radiated power P tot in the case of a broadband emitter as a function of 
the numerical aperture of the collection objective (a) without additional filter (b) with 
a AA = 40 nm band-pass filter. The calculations have been carried out for a single 
HR mirror and for microcavities with three values of N and different values of the 
microcavity mode - emitter wavelength detuning i5. The results are normalized for the 
value corresponding to the HR mirror and the maximal numerical aperture (N.A.=f). 



use larger numerical aperture collection objectives with N.A. close to unity. Note that 
the collection efficiency enhancement is improved by the band-pass filter as shown in 
figure [51(b). This effect comes from the increase of the coherence length due to the 
coupling of emitted light to the microcavity mode, similarly to what has already been 
observed in microcavity LED |54j. We can underline furthermore this phenomenon by 
calculating the enhancement ratio P tot /P where P is the radiated power by the HR 
mirror at the same N.A. than the one used to calculate P tot . The corresponding results 
are shown in figure El (a) for the full spectrum and (b) for the filtered spectrum. We 
note that the maximal collection efficiency is reached for N.A.> 0.6. For N.A. =0.6 and 
5 = 30 nm detuning the maximal enhancement is increased by 1.7 for the full bandwidth 
detection [figure [6j( a)] and by 2.75 for the filtered detection [figure |6j(b)]. 

At low N.A.< 0.6 and filtered detection, the collected power is reduced in the 
microcavity case (N = 4) with respect to the mirror case for detuning larger than 
5 = 30 nm. This situation is reversed if one uses smaller detuning values as can be seen 
for 5 = 30 nm and 5 = in figure EJb). Indeed for low numerical apertures a large 
positive detuning has a detrimental role since the cavity resonance would be obtained 
for too large angles of incidence, out of the collection aperture. 

Finally we have also calculated that for N.A.< 0.82 and the configuration N = 4 / 
5 = 30 nm the total collected power in the case of the filtered detection is larger than 
the full bandwidth detection one with only the HR mirror. 

Note that the evolution of P tot vs N.A. obtained for the optimal NV color center 
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Figure 6. Enhancement factor of the total collected power P to t vs. the N.A. in the 
case of the microcavity with respect to the HR mirror configuration, for which the 
collected power is P (a) without additional filter (b) with a AA = 40 nm band-pass 
filter. The power P to t is calculated in the same way than in figure O 



dipoles orientation holds for all possible orientations. For broadband emission, in the 
worst configuration O = 0, and in the case of a numerical aperture N.A. = 0.6, the same 
Ptot calculation methods yields a total collected power equal to 53% of the one of the 
best configurationin the case of the HR Bragg mirror. Similar results are obtained in 
the case of a N = 4 microcavity, for which the total collected power in the worst dipole 
orientation case corresponds to 50% of total power in the optimal case. 

As a conclusion the simulation shows that the collection of photoluminescence from 
a single NV~ color center can be increased by its coupling to one low order mode of 
a planar microcavity. Even in the case of broadband emission the use of a moderate 
finesse (J 7 ~ 60) can enhance the collection efficiency. In addition such coupling also 
lead to a narrowing of the emitted spectrum and thus to an increase in the coherence 
length of the emitted light. 

3. Experimental measurement of the spectral properties of the NV 
color-center coupled to a single mode of the planar microcavity 

3.1. Experimental setup 

The single emitter studied in this work is a NV - color center in type lb (natural 
specified content in nitrogen 200 ppm) diamond nanocrystals (Micron+ MDA 0- 
0.1 from De Beers/Element6) with a size smaller than 100 nm, obtain by a high- 
pressure and high temperature synthesis. In order to increase the natural content 
of NV center in the nanodiamonds, the particle were irradiated by a high-energy 
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Figure 7. Scheme of the setup used to study the fluorescence from a single NV color 
center in diamond. Oi: microscope objective (long working distance, x40; N. A. =0.60); 
O2: microscope objective (xl6; N. A. =0.32); PH: pinhole of 50 /Ltm of diameter; BS: 
beam splitter plate; Fl: glass filter (Schott OG550) transmitting A > 550 nm; F2: 
interference filter transmitting A > 580 nm and removing the remaining light at 
the pump wavelength; BSC: nonpolarizing 50/50 beam-splitter cube; APDs: Silicon 
avalanche photodiodes in photon counting regime; IRF: interference filter suppressing 
near infrared light at A > 740 nm, which is combined with two GFs: glass filters 
(Schott KG4), resulting in the suppression of the optical crosstalk between the APDs 
due to their breakdown flash [55]. TAC: Time to Amplitude Converter connected 
to a multichannel analyser board allowing to extract the time intensity second order 
correlation function g^ 2 '. 

electron beam to generate vacancies and then annealed during 2 hours at 800 °C under 
vacuum. Then, they were cleaned by oxydative chemicals and dispersed in a propanol 
solution containing a polymer at low concentration (polyvinylpyrrolidone, 1% w/w). 
The mean hydrodynamical size of the nanodiamond in the colloidal suspension is 90 nm 
as measured by dynamic light scattering. The nanodiamond-polymer solution is finally 
spincoated onto the HR dielectric mirror forming a 25 nm thick polymer layer with a 
surface concentration of about one nanocrystal per 10 /im 2 . The microcavity is obtained 
by placing a LR mirror (N = 4 reflectivity, Rlr ~ 90 %) in front of the HR mirror 
(reflectivity Rhb ~ 99.99 %) at a distance that can be adjusted using micrometric screws 
and piezoelectric actuators. From an experimental of view, the actual values of layer 
thicknesses for the HR mirror are: h 2 = 75.3 nm and I2 = 117.1 nm. 

Figure [7] is a schematic of the confocal optical scanning microscope [56] used to 
study the fluorescence of a single NV center in a dielectric microcavity. The light at 
514.5 nm wavelength from a cw argon-ion laser is focused on the nanodiamond-polymer 
layer through the LR mirror, with a long working distance microscope objective Oi 
having a N.A.=0.6. The objective Oi is corrected for the fused silica substrate of the 
LR mirror. Nevertheless, its optical characteristics can be switched, which allows it to 
be used without the additional LR mirror. The LR mirror has a reflectivity less than 
20% in the range 500 — 550 nm allowing for efficient excitation at the laser wavelength. 
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The sample is raster scanned with a piezoelectric stage. The fluorescence from NV 
centers is collected through the same objective Oi, the detection volume being defined 
by a confocal pinhole. The collected light passes through a dichroic mirror and long-pass 
filters are used to remove any residual excitation light. The fluorescence light intensity 
is equally divided into a first beam for spectral analysis by an imaging spectrograph 
(composed of a grating and a cooled CCD array for spectrum record) and a second 
beam for measuring the photon statistics using the standard Hanbury-Brown and Twiss 
(HBT) coincidence setup. The latter measurement allows us to identify the emission of 
a single color center yielding a perfect photon antibunching. The white light source and 
the microscope objective O2 are used to measure the microcavity transmission spectrum. 
Note that in our experimental configuration we obtain a transmission spectrum T(A) 
averaged over the numerical aperture of O2 and over the two polarizations TE and TM 



as detailed in |Appendix A 



3.2. Narrowing of the spectrum of a single NV color center coupled to one 
microcavity mode 

To identify well isolated fluorescent emitters, we first raster scan the sample. For each 
fluorescent spot, the uniqueness of the emitter is then checked by the observation of 
antibunching in the photon statistics of the fluorescent light. Since after the emission of 
a first photon, it takes a finite time for a single emitter to be excited again and then it 
spontaneously emits a second photon, the antibunching effect appears as a dip around 
zero delay r = in the normalized intensity I(t) time autocorrelation function: 

9 «(r) S W)/(t+ 2 T)) . (6) 

The HBT detection setup allows one to record the histogram of time delays between 
two consecutively detected photons on each APD. Considering our detection efficiency 
of a few percentages, this histogram is a very good approximation of the intensity 
autocorrelation function g( 2 >(r). 

Figure M shows different emission spectra measured for the same single NV center 
while changing the cavity length. The first spectrum at the bottom is obtained for a 
closed cavity (i.e., HR and LR mirror into optical contact). In this configuration there 
is no resonant mode in the Bragg mirror stop-band and the fluorescence of the center is 
inhibited. 

By opening the cavity we increase the resonant wavelength corresponding to m c = 3 
and we obtain a resonant mode whose wavelength is located inside the emission spectrum 
of the NV - center. We observe that the emitted fluorescence intensity is maximal with 
a noticeable shrinking of the spectrum, when the resonant wavelength Ao is slightly 
detuned from the free space NV - center emission maximum A^y x (see third curve 
starting from top of figure ED- 

Figure [9j (a) shows the detailed analysis of one of the optimal configurations. The 
experimental transmission T(A) has been measured and this optimal case corresponds 
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Figure 8. Emission spectra of a single NV - center for different values of the air gap 
do between the LR and HR mirrors. From the bottom of the graph to the top the 
cavity length is increased using piezoelectric actuators. We start from a configuration 
where the two mirrors are almost in contact, and pass trough a distance bewteen 
the mirror yielding a maximum intensity with a narrowing of the emission spectrum 
compared to free space emission. The small peak at 607 nm is the two-phonons Raman 
scattering line of the diamond nanocrystal matrix as observed in previous work [57] . 
corresponding to a Raman shift of 1480 cm -1 , different to the one of 1335 cm -1 in the 
case of bulk diamond crystal, probably due to the nanocrystalline form. 



to the situation where the cavity resonant wavelength is almost tuned at the emission 
maximum. Here, the optimal detuning 5 is almost zero since the experimental value of 
N.A. is only 0.6. Figure [91(b) displays the cavity emission spectrum and the HR mirror 
emission spectrum (used as reference) recorded in the same experimental conditions. 
In particular the saturation of the single NV center has been reached for the two 
configurations. Finally figures [9j(c) and[9j(d) are time intensity second order correlation 
functions associated to the emission of single NV center, for the mirror and cavity 
configurations respectively. The vanishing dip at zero delay r = in the g^ function 
insures that we address single color center in both cases. 

The fit of the transmission spectrum T(A) of the cavity as been obtained only 
adjusting the value of the air gap thickness (here d u = 124.1 nm and thus 5 = 1.2 nm) 
governing mainly the resonant wavelength value. As mentioned in Appendix A we 



consider that the numerical aperture value is limited by the illumination objective 
(N.A. =0.32). The spectral width of the calculated linear transmission spectrum T(A) is 
in good agreement with the experimental results. For the comparison between the 
experimental and theoretical fluorescence spectra, we had to take into account the 
characteristics of our confocal microscope. Indeed, it is not possible to reproduce the 
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Figure 9. (a) Transmission spectrum of the microcavity obtained with its illumination 
by a white light source at N. A. =0.32 (see figure [7]). We also show a calculated 
transmission spectrum obtained for do = 124.1 nm. (b) Emission spectra of a single 
NV centers on a HR mirror (grey plot) or inside a cavity (pink plot). The calculations 
are normalized using the microcavity spectrum. We give two sets of calculated spectra 
i) dash curve: without taking into account the pinhole (PH) spatial filtering ii) full 
line: taking into account the confocal PH filtering and a negative defocusing. The small 
peak around 607 nm of graph b) is a Raman scattering line of the diamond matrix. 
Time intensity correlation function displaying photon antibunching in two cases: 
(c) single NV center on a HR mirror (d) single NV center inside the microcavity. The 
fits (blue lines) of the correlation functions have been carried out using the model given 
in reference 1581 . 



width and the amplitude of the experimental spectra only using the expression of the 
overall radiated power dP tot /d\ derived from equation (E]). This is illustrated by the 
dash theoretical curves of figure El (b). We have normalized the theoretical calculations 
using the maximum of the microcavity spectrum. In this case the theoretical HR mirror 
emission spectrum, normalized in the same way, does not fit the experimental data. The 
confocal aperture spatial filtering must be taken into account to obtain a good agreement 
between the theory and the experimental results [59J. The model used to simulate the 



propagation across the microscope is described in Appendix B After fixing the value of 
g?o obtained from the transmission spectrum fit, the defocusing 5z (see Appendix B) is 
the only adjustable experimental parameter. Note that we have already mentioned that 
we have experimentally selected the brightest emitters whose the two orthogonal dipoles 
are parallel to the mirror plane. Consequently we assume in the following calculations 
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Table 1. Counting rates of a single NV color center in the case of the two studied 
configurations: HR mirror alone and microcavity (HR mirror + LR mirror on top). 
Note that the measurements were done on two different single emitters. 



Configuration 


Excitation laser 


Total photon 


Counting rate spectral 




power 


counting rate 


density around Ao 


HR mirror 


11.5 mW 


140 kcounts/s 


0.5 counts/s/nm 


Microcavity 


10.2 mW 


63 kcounts/s 


1.0 counts/s/nm 



that = 7r/2 and a = as it was the case in the theoretical part. The diameter of the 
confocal aperture D = 50 /iin is well adapted for the case of an emitter deposited on 
the HR mirror. In this case, the pinhole does not truncate the electric field distribution 
and we obtain an almost perfect collection efficiency for 5z = 0. Calculations show that 
for an emitter inside the cavity, the pinhole diameter is not perfectly adapted. The 
collection efficiency is thus reduced with respect to the case of the HR mirror. This 
detrimental effect could be reduced by using a larger confocal aperture. The effect of 
the LR mirror on the imaging properties of the system is taken into account by using 
a negative defocusing [61)] • The theoretical calculations given in figure [9j (b) have been 
carried out using a defocusing 5z = —1500 nm almost corresponding to the optimal 
value for the overall signal collection. In this case, we obtain a very good agreement 
simultaneously for the shape of the spectrum and for the relative amplitude between 
the reference (HR mirror) and the microcavity. 

We show that the spectrum of a single NV center can be tailored using a planar 
microcavity in a similar way than for single molecules [38]. From a qualitative point of 
view, the fluorescence spectrum is asymmetric and shows profile steepening from its red 
side which is in good agreement with other experimental observations in microcavity 
LEDs [61]. We also show a good agreement between the experimental results and the 
theoretical calculations taking into account the optical and geometrical properties of 
our experimental setup. 

We then compared the performances of the sole HR mirror single-photon source, to 
the one with the emitter inside the microcavity. Tabled] reports the total counting rates 
(sum of both APDs), and the counting rate spectral density for the two configurations 
at similar input excitation laser power. We measured 140 kcounts/s for the HR mirror 
configuration and 63 kcounts/s in the case of the microcavity. We thus observe a 
reduction of 55% between the HR mirror and the cavity case. The theory gives only a 
reduction of 27% due to the pinhole effect. The difference might come from the filtering 
of the parasitic light (e.g. Raman scaterring) in the case of the cavity. Note that we 
would have obtained an augmentation of 70% without taking into account the pinhole 
filtering. We however observed an interesting modification of the fluorescence in the 
cavity configuration: due to the narrowing of the emitted spectrum from B m 75 nm to 
B Pd 22 nm, the counting rate spectral density is twice larger in the microcavity case, 
which proves that the emission modes of the nanocrystal are to some extent coupled to 
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some microcavity modes for which the fluorescence rate is enhanced. 
4. Conclusion 

We have shown theoretically using a simple dipole model embedded in a multilayer 
system that the fluorescence collection of a single NV diamond color center can be 
enhanced by a factor of two when it is placed into a moderate finesse (J 7 ~ 60) 
microcavity with respect to its emission when it is located onto a single HR mirror. 
The calculations have been done in practical configurations of either Quantum Key 
Distribution corresponding to a room temperature operation and thus broadband 
emission of the NV center, or of Quantum Information Processing associated to cryogenic 
temperatures and monochromatic emission at the zero-phonon line wavelength. We 
also performed room temperature experiments using NV color center in diamond 
nanocrystals and studied the effect of a planar microcavity on the detected fluorescence 
light from a single NV center. We have shown that the fluorescence spectrum can be 
narrowed using such a planar microcavity and that the counting rate spectral density 
can be increased compared to the use of a single HR mirror. Such simultaneous spectrum 
narrowing and fluorescence spectral density enhancement brought by the microcavity 
represent a significant increase in the performances of the diamond NV color center 
based single-photon source in the prospects of open-air QKD applications. 
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Appendix A. Linear transmission calculation in the case of a focused 
illumination 

We give here the expression of the averaged transmission spectrum T(A) of the 
microcavity used to reproduce the experimental results. We used the standard linear 
transmission T(A, 9, P) calculated for a plane wave propagating in the 6 direction with 
a polarization P = {TE, TM} using the transfer matrix method to evaluate the 
quantity: 

T ( A ) = ^ Ff / t T ( A ' TE ) + T ( A ' ™ 1 sin m - A -l) 

2{L — COS um ) Jo 

We assume that the beam limitation only comes from the numerical aperture of the 
microscope objective O2 used to illuminate the microcavity with the white light as 
shown in figure [7] and thus sin 6*m = 0.32. 
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Appendix B. Theoretical description of the propagation in the imaging 
system 



In this section we give the method used to take into account the effect of the confocal 
aperture (diameter D) and the defocusing on the collection efficiency. The defocusing 
is taken into account by adding a perturbation to a reference wave surface propagating 
in our stigmatic optical system as represented in figure IB1I [62| [63] . We consider here 
three wave surfaces: i) £/ which is the object focal sphere of radius / ii) a plane wave 
surface E located in the intermediary space between the microscope objective (focal /) 
and its tube lens of focal /' and iii) the image focal sphere E// of radius /'. We assume 



NV 
center 




PH 

plane 



Microscope 
optics 



Figure Bl. The reference focal sphere X/ is centered on the emitter. This surface 
is transformed into a plane £ between the objective and the tube lens. In presence 
of defocusing Sz', the reference sphere is changed by adding a phase term ipdd- The 
confocal aperture (PH) of diameter D is located at the center of the unperturbed 
reference focal sphere £ /' . Both refractive indices of object and image spaces are 
equal to 1. 



a perfect imaging system and these three surface waves are stigmatic. First we calculate 
the field Ej(q) (where we note q = k + /fc) in the 9 direction on the reference sphere Ej 
using equations (OB). Then we obtain the electric field on the plane wave surface E [64J: 



EfefoO = —L^E/fa. = v/f,Qy = (B.l) 
vcos u 

where rj and £ are the coordinates of the intersection between the emerging ray and E. 
For a given direction characterized by 6' and the unitary vector q', we can write the 
emerging field on the surface E // : 



E f (q') = ^-E s (r, = = q' y f)e^ e '\ (B.2) 

The additional phase term (yjdef accounts for the image space defocusing 8z'\ 

Pdrf(0') = ^(l-cos0')- (B.3) 
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For a given object aperture angle 6, the corresponding image angle 9' is calculated using: 
p = f sin 6 = f sin 6' where p = a/ 7 ? 2 + £ 2 [64J • The defocusing value in the image space 
5z' is connected to the defocusing in the object space 5z using the sine and Herschel 
conditions noting that /'// = Gt where Gt = 40 is the magnification of the objective 
Oi. In the image space, the numerical aperture is weak and the paraxial approximation 
can be used. In particular, the electric field E in the focal plane (x', y') can be calculated 
using a Fourier Transform (FT): 

Ep H (x7 A, y'/X) = FT [E f (rj/f, £//')] . (B.4) 

We can thus obtain the intensity distribution in the image plane by keeping only 
the x' and y' components of Eph [601 E2] which allows the collected spectral density 
power to be calculated by an integration over the pinhole area: 



dPtot , , \ ceo 



dX y ' 2 




PH 



(A) = ^// /(A) \E Pn , x ,(x'/X,y'/X)\ 2 (B.5) 



+ \EpB., V '(x'/X,y'/X)\ dx'dy'. 
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